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We discuss a feedback mechanism between electronic states in a double quantum dot and the 
underlying nuclear spin bath. We analyze two pumping cycles for which this feedback provides a 
force for the Overhauser fields of the two dots to either equilibrate or diverge. Which of these effects 
is favored depends on the g-factor and Overhauser coupling constant A of the material. The strength 
of the effect increases with A/Vx, where Vx is the exchange matrix element, and also increases as 
the external magnetic field Bsxt decreases. 
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Hyperfine interaction with the host nuclei in nanoscale 
GaAs systems, while relatively weak, can nevertheless 
limit the electron coherence time and thereby compli- 
cate strategies to implement quantum information and 
quantum computing schemes in these systems [TJ |2l [31 14] . 
Conversely, ever-increasing control of angular momentum 
transfer between electrons and nuclei in a range of mate- 
rials enables numerous applications precisely because of 
the environmental isolation of the nuclear system. These 
include applications to quantum information processing 
employing NMR [5j . From the perspective of fundamen- 
tal physics, experiments on few-electron systems with 
controllable coupling to the nuclear many-body system 
uncover a fascinating arena of new phenomena with ram- 
ifications for theoretical physics and engineering [5]. 

Experiments on double quantum dots with electron 
number N = 2 have uncovered and exploited an intrigu- 
ing phenomenon called the "Pauli blockade" [T in which 
two electrons with parallel spins are forbidden from com- 
bining in one dot by the exclusion principal. In transport 
or in gate pulsing, even when such a transition becomes 
energetically favorable, it can only proceed via a spin 
"flip-flop" process in which angular momentum is ex- 
changed with the local nuclei. Repeating the spin trans- 
fer modifies the character of the nuclear spin distribution. 
One metric for the nuclear state is the difference between 
the total Overhauser fields of each dot. These are the ef- 
fective Zeeman fields which the electrons experience due 
to nuclear polarization. Several recent experiments ad- 
dressed transfer of angular momentum from the electron 
system to the nuclear bath through various pumping cy- 
cles. One experiment claims that under a specific, re- 
peated pulsing sequence (see below) [8] the polarizations 
in the two dots tend to equilibrate; a phenomenon which 
has been numerically reproduced [SJ However, an- 
other similar experiment claims to find a large difference 
induced between the Overhauser fields of the two dots 
[TT]. The theory which we describe here does not claim 
to explain either experiment. 
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FIG. 1: Top: electronic states near the (1,1) to (0,2) stabil- 
ity diagram transition. Bottom: overlap of the ^ state with 
5(0,2) (a), with |L|R|) (b) and with |L|Rt) (c). Parame- 
ters: Bext = 0.2 r, 7 = 1.2iieV, A = 1000, Ec = O.GmeV, 
Vx = IfieV. 



Here we describe a force toward either equalizing or 
inducing differences between the Overhauser fields in the 
two dots. The direction of this force depends on the 
spin of the initial electron state (i.e. the direction of the 
electron-nuclear "flip-flop" process) as well as on the sign 
of the product of the g-factor g of the host material and 
the sign of the Overhauser magnetic field of the nuclei 
of Ga and As as compared to the direction of the nu- 
clear spin (they are anti-parallel). Assuming GaAs, we 
describe two pulse sequences which differ in the choice of 
the initial electron state, which consequently have a force 
tending to cause the Overhauser fields in the two dots to 
equilibrate or to diverge. 

Electronic States of the Double Dot with N—2 - We cal- 
culate the electronic states of the two electron {N = 2) 
double dot within the Hund-MuUiken formalism [T5] de- 
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veloped for the hydrogen molecule. Typically, in this 
method, eigenstates of total spin, singlets and triplets, 
are employed as basis states. However, since we wish 
to study the inhomogeneous Overhauser effect due to 
different effective magnetic fields in the two dots, we 
choose a basis which diagonalizes, at the single particle 
level, the z-component of this inhomogeneous field and in 
which the spatial dependence of nuclear spin flips induced 
by electronic spin "flops" is transparent. The basis is: 
Un} = {IRtRi) , ILtRx) , IL^Rt) , iLtRt)}, where L and 
R indicate the orbital states of the left and right dot, the 
arrows denote spin direction |13j . Two remaining states 
of the Hund-MuUiken model, jL^Rj^) and |L|L|) are not 
relevant to our analysis. Note that |R|R|) is the stan- 
dard S{0,2) state and |L|R|) is the standard |T+) state. 
The hyperfine Hamiltonian for two electrons is properly 



written: 



Hhf = -^^[(5(ri-R,„)SrI™®l+l(»<5(r2-R„)S2-I,„] 

m 

(1) 

where and Si are operators in the subspace of electron 
i (first quantized representation) and m is summed over 
a total of M nuclei (typically M ~ 10^); and where v 
is the volume per nucleus. We assume, for simplicity, a 
single nuclear species with spin f/2. Then, constraining 
the maximum Overhauser field to be 5.3 T [T3j leads to 
a coupling constant A = — 270 fieV. We incorporate the 
matrix elements of Hf^f from Eq. [l]in our basis {Cn} into 
the Hund-Mulliken Hamiltonian which gives 
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which is correct to leading order in {L\B), where e is the 
potential "detuning" (the difference between the electro- 
static potential bottom of the left and right dots) and 
where we have taken the orbital energies of L and R to be 
zero for simplicity. We include only two Coulomb terms: 
the charging energy Ec = Vrrrr - Vrlrl and the ex- 
change matrix element Vx = Vlrrl ilSj . Also in Eq. [2] 7 
is the tunneling coefficient; Ez = gfJ-sBext is the Zeeman 
energy for a presumed external field Best, which defines 
the z-axis of the problem, with fis the Bohr magneton. 
The lower left hand side of the matrix is the complex con- 
jugate of the upper right hand side. Note that the matrix 
elements of H in this electronic basis remain operators 
in the Hilbert space of the nuclear coordinates [T3]: 

A 

^''^ ^""m^ C(Rr„)V'/3(R™)4 (3) 

m— 1 

where a, (3 S Note that previous researchers 

have typically ignored the transition term I^", which 
we see from Eq. [2] can lead to a direct transition between 
|R|R|) and |L|R|) and causes a spin flip in the barrier. 
This term, and the other overlap terms (e.g. cx {L\R)), 
could be significant in the case of large Bf,xt and small 
7, i.e. where |R|Rx) and |L|R|) anti-cross deep in the 
(0,2) regime. However we will henceforth ignore hyper- 
fine terms in Eq. l2]proportional to wavefunction overlap 
(e.g. {L\R) and ) and leave exploration of the barrier 
nuclear spin flip to a later publication |17j . 

Nuclear spin flip location - The crucial feature of Eq. 



|2]is that the |L|R|) state is coupled to |L|R|) via a term 
which flips a nuclear spin in the right dot [if") and it is 
coupled to |L|R|) by a term that flips a nuclear spin in 
the left dot {li^^) ■ In the absence of flip-flop coupling to 
the |L|R|) state, the upper left 3x3 matrix in Eq. |2](see 
also yellow highlighted region of Fig. 2) has a ground 
state, which we denote: 

\^) - a{e) IRtRx) + h{e) iL^Rx) + c(e) IL^Rt) . (4) 

As shown in figure 1, at large (positive) e, Y^) 
|R|R|) = S{0,2) and at large negative e, |^) becomes 
an unequal superposition of |L|R|) and |L|R|). Even 
when Vx > \{lf" — lf^)\, the inhomogeneous Overhauser 
effect will produce a preference for either the |L|R|) or 
the |L|R|) component of \E' (see figure 1), with the elec- 
tron down spin preferentially located on the dot with 
smaller J^. In the first electron pulsing sequence which 
we describe, the electron state is initialized at large e into 
1^) « S{0, 2) and detuning is swept approximately adia- 
batically through the ^' - |L|R|) anti-crossing. The posi- 
tion of this anti-crossing, e, is determined by the energy 
of |L|R|> EE |T+), denoted E{T+) (see figure 1), which is 
determined by Bext- Insofar as b(e) 7^ c(e), a transition 
from ^> to |L|R|) will preferentially induce a nuclear spin 
flip (down) on the side with the larger 7^. This tends to 
equilibrate the values of if" and if'". In the second 
pulse sequence the electrons are initialized into |L|R|) 
and the state then transitions to The same feedback 
mechanism preferentially now causes nuclear spins to flip 
up, but still on the side with the larger Iz, thus leading 
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to a tendency for — to grow. Both of these 

sequences can be experimentally implemented [HI [18] . 
Our further analysis focuses mainly on the first pulse se- 
quence. 

Note that the preceding argument depends on the sign 
of A which in turn depends on the sign of g and the sign 
of the effective Overhauser field which, for Ga and As, 
are anti-parallel to the nuclear spins ^19^. 

Nuclear States -To further analyze the Hamiltonian, 
Eq. [2] it is helpful to introduce a simplified basis for 
the nuclear states in which all of the nuclei are either in 
the left or right dot and all within a given dot interact 
equally with the electron. In other words, |t/'L(r)p is 
taken as a constant within a spherical "box" of some 
volume, V. In this model, which we refer to as the "box 
model," the squares of the total angular momenta 
are conserved, where = (vA/V) X^mea where 
a G {L,R}. Thus, the electrons essentially interact with 
two composite nuclear spins, one on the left and one on 
the right. The nuclear state basis is {II, Ir, Ilz, Irz} 
(where -|- 1) is the eigenvalue of and laz 

is the eigenvalue of Finally, for given it is 

convenient to transform to the basis of A = I^^ — Irz 
and s = + Irz ■ In this basis the z-components of the 
nuclear operators have non-zero matrix elements on the 
diagonal blocks, but the raising and lowering operators 
connect different (A,s) subspaces (see Fig. 2). 

The strength of the narrowing force depends on the 
ratio r = c/b at e. This depends on A and on B^xt- 
For example, smaller B^xt results in smaller (or more 
negative) e, where, as shown in Fig. 1, the ratio c/h 
increases (for A > 0). Exactly how large c/h can get 
depends on Vx which, in the example of Fig. 1, we have 
set to 1 ^leW |2Q]. 

In Fig. 3 we plot the value of r(e) as a function of 
B(,xt for various values of A. The key point is that r{e) 
increases monotonically with A (cf. yellow highlighted 
region of Fig. 2), however it also decreases monotoni- 
cally with B^xt (and hence e). Interestingly, because the 
state is coupled equally to |L|R|) and |L|R|), 
the value of h/c is independent ofj. 

We note that the flip-flop process naturally also de- 
pends on the rate at which e is swept since, in order to 
be adiabatic and remain on the lower branch of the ^' - 
|L|R|) anti-crossing the e variation must be sufficiently 
slow. More generally, the character of the state evolution 
can be examined as a Landau-Zener tunneling problem 
[2T] or else evaluated numerically [TZ] . 

The evolution of the full nuclear state is complex and 
the experimental manifestations of that evolution are am- 
biguous. Nevertheless, as a possible baseline for more 
detailed studies of the nuclear evolution, we describe a 
simple, incoherent model which results in narrowing of 
the distribution of A. 

If we assume that the system is in the well-defined 
state 1^) (g) {Il,Ir,Ilz,Irz} and the detuning is moved 
quickly to e and held there for time r, we can compute, 
by Fermi's golden rule, the probability for a nuclear spin 
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FIG. 2: Hamiltonian for three sectors of the nuclear difference 
quantum number (A — 1, s — 1), (A, s), (A -|- 1, s — 1). In the 
above, ±A ± 1 and s ± 1 are shorthand for {vA/V){±A ± 1) 
and {vA/V){s ± 1) respectively. 



to flip in the right dot as: 
'^r{IlzJrz IlzJrz - 1) = rfl(s, A 



s-l,A-l) 



^J^\{lLzlRz-l\l'y'\lLzlRz] 



(5) 



where we have suppressed the II , Ir dependence for 
brevity and where the matrix elements of the ladder op- 
erators are given by the well-known formulas: il.a± = 

± l|/± = y/o(/a + l)- /c.z(/a.±l). 

Similarly, the flip probability in the left dot is propor- 
tional to the c component of ^' 



ri(s,A^s-l,A + l) 



(6) 



If we denote the probability distribution for the nuclear 
state (at fixed Il,Ir) as W{s,A), then the condition for 
W to be stable in its dependence on A can be written 
(cf. Fig. 4a): 



W{s + 1, A + 1)Tl{s + 1, a + 1) = W{s, A)Tr{s, A) 

W{s,A + l) = Wis,A)^,^ (.,A + l)|c(A + l)P 

(7) 

where we have assumed that W{s) « W{s + 1) and we 
have used the fact that b and c depend very weakly on s 
(only through the s-dependence of e). 
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butions of II and Ir, and hence will show no narrowing 
of VF(A). Thus our box model can qualitatively explain 
the narrowing effect or the saturation, but not both. 

We believe that a full understanding of these phenom- 
ena depends on the variable coupling of the electron wave 
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FIG. 3: The wave function ratio r = c/b evaluated at the 
^ — r+ crossing point, e, as a function of B^^t for various 
values of the nuclear spin z-component difference A. r(e) is 
monotonically increasing with A and decreasing with Bext- 



Recursion relation Eq. [7] can be solved iteratively 
and the influence of the narrowing force evaluated. In 
Fig. 4 we have plotted W(A) computed with the ratio 
riij-Zr^L- set to unity to show only the narrowing from 
the inhomogeneous Overhauser effect described here with 
the same electronic parameters as in Fig. 1, and with 
II = Ir — 1000; including the fi's induces more narrow- 
ing. For comparison we show the T — > oo thermal distri- 
bution of A, averaged over s, also for Ii^ — Ij^ — 1000. 
Inset (a) shows the ratio of the root-mean-square (rms) 
A in the thermal distribution, ctt, to the rms A with 
the narrowing force at varying Bext, cr{Bext)- A substan- 
tial narrowing of VK(A) results from the inhomogeneous 
Overhauser effect. 

Discussion - Experimentally, the \^!) to |L|R|) pulse 
sequence polarizes only about 1% of the nuclei, even 
when running sufficient cycles to flip all of the nuclei 
|18) . This saturation of the nuclear polarization is still 
an open problem. A recent article by Yao discusses a 
model similar to that described herein. In that paper, no 
mechanism for stopping the flip-flop process is proposed 
when the pumping continues (as it does in experiments) 
beyond ~ 10^ cycles. In our model, polarization will sat- 
urate when both 1^^ = —I^ and Ifj^ = —Ir, implying 
that s = —II — Ir- However, the resulting distribution 
of A will then mirror the difference in the initial distri- 
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FIG. 4: (main) Reduced distribution W{A,s), calculated 
from Eq. [t] (solid lines), for 14 = 1 /ieV as a function of A 
for various B^xt = 0.05,0.10, ...,0.75T (lower fields have nar- 
rower W); and thermal VF(A) (dashed), averaged over s, all 
with II = Ir = 10^. Inset (a) narrowing factor ctt /o-{Bext) 
versus B^xt- Inset (b) Illustration of Ilz — Irz plane. A and 
s are the diagonal coordinates, with A = Ilz — Irz- 



function to different groups of nuclei, so that conserva- 
tion of the magnitudes of two spins, II and Ir, is not 
required. Such a model with multiple interacting com- 
posite nuclear spins, incorporating the narrowing effect 
described here as well as the Landau-Zener tunneling be- 
havior near e, in some parameter regimes shows the po- 
tential to send |A| — s- while reducing the spin flip prob- 
ability, slowing the growth of total polarization; for other 
parameters, |A| grows large despite the narrowing force 
described here [21]. 
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